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>> \ , We obtain the photon spectrum induced by a cosmic background of unstable neutrinos. We 

study the spectrum in a variety of cosmological scenarios and also we allow for the neutrinos having 
a momentum distribution (only a critical matter dominated universe and neutrinos at rest have been 
considered until now). Our results can be helpful when extracting bounds on neutrino electric and 
magnetic moments from cosmic photon background observations. 
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' I. INTRODUCTION 

Recent indications for neutrino masses |l| have strengthened the interest in physical effects linked to massive 
neutrinos. In this paper we shall be concerned with one of these effects, namely, with the fact that massive neutrinos 
Q\ | can have non-zero electromagnetic form factors and can be unstable due to radiative decay processes. 

Specifically, we shall be interested in the contribution of decays v% — ► fjj, when coming from a cosmological Vi 
density, to the cosmic photon background. A previous study of such phenomena was performed by Ressell and Turner 
Q (see also ||). These authors considered the photon background at different photon energies and constrained the 
radiative decay of a cosmic density of massive neutrinos. Recently, Biller et al. have improved substantially the 
constraints in the infrared region by making use of the effect of 7 — 7 interactions giving electron-positron pairs, 
where one 7 is a background infrared photon, and the other is a TeV 7 coming from an active galaxy. We may expect 
further improvements in the future, along the lines used in Q or perhaps with some new ideas and/or observations. 
Constraints on the neutrino lifetime lead to limits on the magnetic and electric transition moments, fiij and ey, of 
. £h neutrinos. The neutrino lifetime r is related to these parameters, 
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where Am 2 = mf — mj, and B is the branching ratio BR(vi — > 1^7). 

All the studies performed until now make the simplifying assumptions that the decays are produced 1) in a critical 
matter dominated universe and 2) by neutrinos at rest. Although this might be enough for the results obtained until 
now, in the future it may be interesting to study the contribution to the photon background with more generality. In a 
more general framework we can be able to see how the limits on the neutrino lifetime depend on different assumptions, 
and we can extract more precise and reliable bounds. Also, in the eventual case of a positive signal that may come 
from decays of relic neutrinos, the general study of the consequences of the decay that we perform in the present paper 
would become quite necessary. Apart from the above mentioned references, other papers treating unstable neutrinos 
in a cosmological context are listed in ||. 

Non-vanishing values of fiij and/or lead to other potential effects besides the ones considered in this paper. 
For example, plasmon decay into neutrinos in stellar media is constrained by arguments of stellar energy loss. This 
leads, in general, to limits that are quite stringent [||. We would like to point out that when neutrinos are close to 
a degeneration in mass, the limits obtained using Eq. (p]) may be relevant since in this case Am 2 /^ < 1. In any 
case, the spirit of the present paper is not phenomenological. Rather, as we said, we would like to generalize previous 
studies on the subject. 
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The organization of the paper is as follows. In Sec. y we examine the contribution to the photon background of 
the decay Vj — > z/j7 in different cosmological scenarios, keeping the assumption that the neutrinos i>\ decay at rest. 
In Sec. M we drop this assumption and study the decays of neutrinos with a momentum distribution. We devote 
Sec. to the conclusions. Some technical details are developed in two appendices. 



II. DECAYING NEUTRINOS AT REST 



In this section we shall calculate the photon spectrum produced by the decay of a cosmic background of neutrinos 
at rest. We shall consider the two-body decay of a neutrino with mass m, into a photon and a neutrino with mass 
rrij, with wn > mf vi Uj + 7. The subscripts i,j = 1,2,3 stand for any neutrino mass eigenstate, which is a linear 
combination of the three weak eigenstates v e , v^, v T . When a neutrino decays at rest the photon energy is given by 

Ato 2 

€0 = W (2) 

with Am 2 = mf — to 2 . The cosmic expansion redshifts the photon energy. A photon that at present has energy E 
was produced at a redshift Zq given by 

! + *> = §■ (3) 

Let Fe be the present energy flux of photons with energy E produced by neutrino decay. The flux per unit energy 
and solid angle is given by 

d2p E = E d 2 F n 

dEdQ dEdn' { ' 

where F n is the particle flux at present. It is related to the particle flux at emission time by the relation 

d ~£r = (TT^^ (<(zo)) = h (TT^F^ ( ' (zo)) ' (5) 

where we have included the factor of dilution (1 + zq)~ 3 produced by the expansion of the universe. The photon 
density emitted at z is given by the usual decay law 

6n 7 (t(zo)) = B— n Vi (t(zo)), (6) 

T 

where r is the neutrino lifetime and B is the branching ratio for the radiative decay. For a fixed emission time 
St = H- 1 (z )dE/E, where H(z) is the Hubble expansion rate at time t(z). Writing everything together we obtain 

dEdn 4tt (1 + zo) 3 tH(zo)' [) 
Choosing a time t p <ti t, otherwise arbitrary, and calling the expansion age of the universe to, we can write 

"^(*( z o)) n Vi (t p ) t(z )-t p _ t(z Q )-t p 

JTTW = (i + ^fa)) 8 exp — 55 niH (<o) exp — ' (8) 

where n Vi (to) would be the present number density of neutrinos if they did not decay. Taking t p = and substituting 
Eq. (||) into Eq. (Q) , we finally obtain 

d 2 F E B h Ui (t Q ) t(z ) 

exp , (9) 



dEdVL 4tttH(z 

for any E < £q. For E > eo the photon flux vanishes because photons cannot be produced with energy larger than eo 
by neutrinos decaying at rest. From now on we shall set B = 1. The expansion time t(z) is given in terms of H(z) by 
the following integral 

f t[z) f°° dz' 1 

tlz> = l dt = L TTPmFi' < 10 > 
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FIG. 1. We plot the normalized spectrum (cPFs/dE <Kl)/{n Ui (to)/47r) versus -E/eo when ri/o = 10. The solid line represents 
the cosmological model Oo = 1 and Qa = 0. The dashed line is the model fio = 0.3 and Q\ = 0. The dotted line is Oq = 0.3 
and Qa = 0.7. 



Equation (|9j) is our final expression for the photon energy flux per unit energy and solid angle produced by the 
decay of a cosmic background of neutrinos with negligible velocities. It holds for any isotropic and homogeneous 
universe and for any equation of state for the cosmic fluid. 

Once the expansion rate H{z) is given, one can calculate t{z) with Eq. ( |Tc| ) and then calculate H(zq) and t(z ) 
using Eq. (^). Plugging the result into Eq. (||) one obtains the photon spectrum as a function of the photon energy E. 
As examples we shall consider three different cosmological scenarios: a power-law expansion and two observationally 
more motivated scenarios, an open universe filled only with matter and a flat universe with matter and a nonzero 
cosmological constant A. 

First, let us take a power-law expansion. The Hubble expansion rate and the expansion time are given, respectively, 

by 

H(z)=H (l + z) n , (11) 

Hot(z) = ""TT - T w' (12) 
n (1 + z) n 

with n > and Hq the Hubble constant [j]]. Particular examples include a critical A dominated universe (n = 0), an 
empty open universe (n — 1), a flat matter dominated universe (n = 3/2), and a flat radiation dominated universe 
(n = 2). For a critical matter dominated universe we obtain 

d 2 F E 1 h Vi {t ) fE\ 3/2 2 1 (E\ 3/2 

exp--— — — . (13) 



dE d£l 47T tHq \ eo / 3 tHq \ eo 

The only case studied in the literature is a critical matter dominated universe (flo = 8irGp m /3HQ = 1, where p rn is 
the present matter density) with rrij = |^| In this case and when r 3> -H^ 1 we obtain 

1 h^to) ( E \ 3/2 ^ 



dEdn 4?r tH \m l /2 
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FIG. 2. Same as Fig. | but with tH = 1 



in agreement with previous results. 

For an open universe with matter density in critical density units f^o < 1 one has 



H{z) = H y/Sl k (l + z) 2 + no(l + z) 3 , (15) 

„ v i y/n k + n (i + z) n + V^hTn^TTTj 

H t(z) — — 5777 m = , (Id) 

n k i + z n 3 k /2 y/W + z) 

being fi^ = 1 — il - 

The last example that we solve is a flat universe filled with matter and a nonvanishing cosmological constant A. 
The contents in matter is given by f2o and the contribution of the cosmological constant is fl\ = A/3Hq = 1 — S7o- 
The expansion rate and the expansion time at redshift z are given, respectively, by 

H{z) = H 0V /{l A + n Q (l + z)3, (17) 
H t(z)= 2 ,-^+y"A + n (l + z)3 

We plot in Figs. |l|, ^, and |^ the photon spectrum (||) in the three cases listed above, for three different values of r, 
tHq = 10, 1, and 0.1, respectively. For an open universe and a flat universe with a nonzero cosmological constant we 
choose the values of the density parameters that seem to be observationally favored: Qq = 0.3 and S7a = 0.7 ||. 

For the first model, fio = 1 and fl A = 0, one can see that the spectrum has a maximum at 

o \ 2/3 /x 2/3 



2 TH ») eo = UJ eo = IT^)< (19) 

in the kinematically allowed region, E max < eo, as long as r < to- The value of the normalized spectrum 
at the maximum is (d 2 Fe / dEdft) / (n Vi (to) /47r) = 3/2e = 0.5518 . . . independent of tHq. As tHq decreases the shape 
of the spectrum becomes narrower and peaked to lower values of E. This reflects the fact that when the lifetime 
r is smaller than the age of the universe to, most neutrinos have already decayed, and the energy of the photons 
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produced in the decay has been redshifted to low values of E. The differences in the spectrum for the different 
cosmological models are easy to understand. Let us compare the critical matter dominated model with the flat model 
with nonvanishing A (the open model falls in-between). For E — ► the photons are produced at large Zq, when the 
Hubble expansion factor for the nonvanishing A model is smaller than H(z) in the critical matter dominated model 
by a factor y/1 — Q\, while the term exp— t/r in Eq. (|]) becomes unity for both models. Therefore, the spectrum 
is higher when A > for E small. On the other hand, when E — > eo, the photons have just been produced, at the 
same expansion rate H but with an age for the A model older than the age for the critical matter dominated case. 
Consequently, the term exp— t/r makes the spectrum for the critical matter dominated model higher than for the flat 
nonzero A model in the range of large E. The position of the maximum in the spectrum when one has a nonvanishing 
A is given by 

*- - (H" G+viX < (H 3/3 - <20) 

which is smaller than the value of E max when A = 0. 

We have shown that the differences in the shape of the photon spectrum introduced by the cosmological energy 
density parameters are small. Only a very accurate determination of the photon spectrum by observation could be 
useful to discriminate among different cosmological models, provided that neutrinos turned out to be unstable with a 
priori known tHq ~ 1. On the other hand, if the photon spectrum was well measured and r was unknown, the small 
sensitivity of the spectrum to the energy contents of the universe would allow a fair determination of r. 



III. UNSTABLE NEUTRINOS WITH A DISTRIBUTION OF MOMENTA 



In this section we shall calculate the spectrum of photons produced by a cosmic background of unstable neutrinos 
which are not at rest but decay with a statistical distribution of momenta fi{pi). Besides the expansion and geometric 
effects discussed in Sec. |n|, the photon spectrum will also include features arising from this momentum distribution. 

As always we consider the two-body decay 
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Vi{ei,^i) — > Vj{tj,Pi) + j(e,p), (21) 

where the energy and three-momentum of each particle are indicated. A photon is produced with an energy e at time 
t(z); at present its energy has been redshifted to E with e = (1 + z)E. 

The present photon energy flux per solid angle with energy between E and E + dE can be written as 

d 2 F E 1 f d 3 p 1 , , s 

1 ' -ES&faz), (22) 



dn 4vr J (2tt) 3 (1 + z) 3 

where d 3 p/(2ir) 3 5f~ l (p 1 z) is the photon number density, including both photon helicities, produced by neutrinos 
decaying at redshift between z and z + dz; these photons are produced with energy between e and e + de. The origin 
of the term 1/(1 + z) 3 is the same as when the i/i's decay at rest. The increment 5f 1 is given by the collision term of 
the Boltzmann equation jjj 

[/.(l - /)')(! + /-,) - (i - U)!lH St. (23) 

The functions fj and / 7 are the statistical momentum distributions of the decay products Vj and 7, respectively. The 
averaged square probability amplitude of the decay ( pl[ ) is given by 

tn 3 1 

|.M| 2 = 16^-^-. (24) 

For a fixed redshift z (or fixed e) the time increment is St = H(z)~ 1 dE/E. The inverse process jVj — ► Vi can be 
neglected when m t 3> TOj since at the late times (r ~ Hq 1 ) we are interested in the photon energy is about 1CP 4 
cV ~ 3 K, much smaller than the neutrino mass rrii ~ 0.1 eV. In order to obtain analytic expressions for the final 
spectrum we will also neglect Pauli blocking and stimulated emission, which are expected to produce small corrections. 
With these approximations we can write 

*f b 1 r d 3 p f rffc f d 3 p- 3 ^ )Ui . ){Pi _ p ._ p)]M? ifi (25) 



dEdn An J (27r) 3 2e J {2i:) 3 2t l J {2ir) 3 2e 1 v ' ^ ^ J r " 1 H(z){l + z) 3 ' 

We shall now relate /, at time t{z) with the distribution function /j that the neutrinos i/j would have at present if 
they had not decayed. The evolution of /$ is governed by the Boltzmann equation 

d h _ 1 f *fi f ^ :(27r) 4, (4)(j)i _ p ._ p)|M|2x 



eft 2e, 7 {2i:) 3 2e J J (27r) 3 2e 

[/,(!- / i )(l + / 7 )-(l-/ i )/ J ./ 7 ]. (26) 

Neglecting Pauli blocking, Bose stimulation, and the inverse process, and using Eq. (pij) we get the differential equation 

dfi(Pi(t),t) lmi , , 

-77 = fi(Pi{t),t), (27) 



which can easily be solved to obtain 



t dt' 

fi(pi(t),t) =fi(pi(t p ),t p )exp- I —t^, (28) 



where we choose the arbitrary constant time t p to be t p <C r. Including the (1 + z) 3 that appears in (25) we can 
write 

Si iPi (*) , *) .3 - ^ _ Si {Pi (tp), tp) ,3 , , /■* dt' m l 



(1 + ^) 3 " W (1 + ^fe)) 3 " Vi " 1 4 r Ci (f) 

= / i («;*o)d 8 rexp- f-^, (29) 
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where q = pi(t p )/(l + z{t p )) — pi(t)/(l + z) is the neutrino momentum at to ■ The distribution function fi(q, to) would 
be the neutrino distribution function at present if the iVs were stable. In appendix ^ we derive its expression. 
Making use of Eqs. @, and (§|) we get 

^-(2n) 4 S^(p i -f j -p)x 



d 2 F E 1 f d 3 p f d 3 q 



dEdtt Ait J (27r) 3 2e J (27r) 3 2 ei J (27r) 3 2 ej 

mS 1 1 , [Uf m 

l^--F77Z\fi{Q, t 0) ex P- / 77a' 

Am 2 t H(z) J t t €i(t') 



6(e t - e 3 - e) Wn^-—J~Mt ) exp- / — -L2-, (30) 



with 



Pi = (1 + z)<f, (32) 
e = |p| = (1 + z)£. (33) 

The momentum pj can be integrated out using 8^ 3 \pi — Pj — The lower limit of the integral in the exponential 
can be set to 0, t p = 0. Assuming isotropy for /, and using Eq. ( |33| ) we can finally write the following expression for 
the photon spectrum produced by the decaying neutrinos 

L 1 dtl{z{t))= / ™ tt ) Uz), (34) 



dEdn Jo Jo l + zH(z) 

where I(z{t)) is the energy flux per unit photon energy and solid angle produced between t and t + dt and is given by 



. 1 f d 3 q ~ lrrii T(y,z) 

/(Z = T / 7^3 fiW ex P 

47T J (27r) ci t e l 



x(l + ^ (^-) 2 s((l + z)E ^V). (35) 

The redshift zo has been defined in Eq. ([|). The angular variable is the angle between the momentum of the 
decaying neutrino and the momentum of the emitted photon, and 

y^- (36) 

rrii 

We have defined the "kinematically contracted expansion time" 

n \- f t{z) ^» m - r dzi 1 i m , 

J o ei{z{ t))-J z l + z'H{z>) + 

which reduces to the ordinary expansion time only if y = 0, T(0, z) = t(z). In appendix |b| we calculate analytically 
the function T(y, z) for the three scenarios defined in Sec. Ill 



Equation (34) has to be multiplied by a factor 2 if, as is usually assumed, there is a cosmological population of 
antineutrinos T>i with the same distribution function as the i>i and the CP symmetry is not broken in the decay. 

The physical meaning of the different terms in Eq. ([35]) is straightforward: there is a sum over all possible momenta 
with which the neutrinos ui can decay, each particular momentum weighed by the decay probability per unit time 
T~ 1 (mi/ei)fi exp —T It. The time dilation factor Ci/rrii that appears multiplying r, in the factor before the exponential 
and in the exponent, takes into account the relativistic delay in the decay of moving neutrinos. The term (1 + z) 
accompanying E stems from the cosmic expansion, which redshifts the photon energy. The term ((1 + z)/(l + zq)) 
accounts for the photon phase space. Finally, Dirac's delta enforces energy conservation. Each photon is Doppler 
shifted by the motion of the neutrino that emits it. The energy of the emitted photon in the rest frame of the decaying 
neutrino eo is blueshifted if cos 9 > (emitter Vi and emitted 7 moving in the same direction) or redshifted if cos 9 < 
(yi and 7 moving in opposite directions). 



Equations (34) and (|35|) are the main result of this section, they give the spectrum of photons produced by the 
decay of a cosmic background of neutrinos. They hold for an arbitrary Robertson- Walker geometry and equation of 
state of the universe, and for any isotropic, nondegenerate momentum distribution of the neutrino background. 
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As a consistency check we evaluate Eqs. ( |34| ) and (^|) when fi(q) = (2ir) 3 n Ui (to)d^ 3 \q), i- e -i au the are at res *i 
and we recover Eq. @. 

We can also calculate the total photon flux per solid angle at present time produced by unstable neutrinos with an 
arbitrary fi(q) 

dF n fdE d 2 F E 1 f d 3 q ~ ( %0)\ .... 

ln^J^dEdn = ^J(2^ m {'-^ —)■ (38) 

If r < Hq 1 we find 



fiW = T~ n "t\ t oh ( 39 ) 



dil 4tt J (2tt) 3 j v ^ 4tt 

as one should expect, since in this regime all Vi have already decayed into photons. 

Next we shall consider the example in which fi is given by the momentum distribution of the standard stable 
neutrino background (see appendix |A|) 

fM = — M L (40) 
exp|gl/T + 1 

where T is the present neutrino temperature and we have set the neutrino chemical potential to zero. We shall study 
the limit 3> To. The motion of the nonrelativistic vi will introduce small corrections to the photon spectrum when 
compared to the spectrum produced by neutrinos decaying at rest. In order to calculate these corrections we shall 
need the asymptotic expressions of T(y, z) when y = \q\/m,i — > 0, which are calculated in appendix [B|. In most cases 
we can write 

H T(y,z)=H t(z)-ay n + ... (41) 

where a is a constant and we neglect terms that decrease faster than y n when y — > 0. Plugging Eqs. (^l|) and ( |40| ) 
into Eqs. (^J) and ([35]) we obtain 

1 n n {t ) _ t(zo) {, , (2 - ^)C(n + 3)T(n + 3)a f T \ n 



dE dQ AtttH(zo) CXP t \ + 3({3)tH \mj + ' ' ') ' ^ 

where C(x) is the Riemann zeta function and r(x) is the Euler gamma function. For instance, in a universe with 
ilo ~ 1 one has n = 3/2 and a = n~ 1/2 r(l/4) 2 /6V^F, see Eqs. (||), ( |B22| ), and ( p26|) . Equation @ is then 

_ 1 MW^Jwf 1 + *('-^>OT 2 fjif + ,..|, ,43, 



dTdfi 47rriJ(2 ) r I 48C(3)^ /2 riJ 

with C(9/2) = 1.0547.. ., C(3) = 1.2021 .. ., and T(l/4) = 3.6256... 
It is also possible to obtain the following asymptotic form of T(y, z) 

H Q T(y, z) = H t(z) +0y n lny + ... (44) 

with [3 a constant, which renders 



d 2 F E 1 n„ 4 (t ) / ( (4-i)C(n + 3)r( n + 3)/J/T 



in, 



dEdQ~ 4nTH(z ) eXP t [ 1+ " 3C(3)rF \m J ln T + ' ' j ' (45) 

For example, in a critical radiation dominated universe we have to make use of Eq. (^]) with n — 2 and /3 = 1/2, see 
Eq. Q), 

^ 1 fi^fa) / + 15C(5) /Jby ln m 1 + ; (46) 



dTdft AtttH(z ) * t \ 2C(3)tH q \m 
with C(5) = 1.0369... 
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In all the above examples the leading correction is independent of E and positive. There are more photons at 
present, with a given energy E, compared to the case when the neutrinos decay at rest. The dominant effect is the 
time dilation in the decay of the neutrinos included in T(jj,z)/t. Photons with present energy E are produced at 
1 + z ~ £o/E; there are more neutrinos available at time t(z) if the decay is delayed, and therefore more photons 
with energy E can be produced. There is a second effect of time delay that goes in opposite direction, the number 
of photons produced at t{z) decreases since the decay probability is inversely proportional to the dilation factor, but 
this second effect is order (To/mi) 2 and therefore subleading. 

We notice that, in general, the natural expansion parameter for T(y : z) is y 1 ! 2 or equivalently (To/m^) 1 / 2 in the 
final photon spectrum, if we assume the distribution (|4(j[). High precision determination of the photon spectrum 
would require analytic calculation of higher orders in (To/mi) 1 / 2 . In addition, the small corrections introduced by 
Pauli blocking and stimulated emission can be important when compared with the relativistic corrections we have 
calculated ||. An accurate determination of the photon spectrum including all these statistical effects can only be 
performed numerically, which goes beyond the scope of the present paper. We think that the analytic expressions 
that we have found can be useful to understand the physics of a late decaying cosmic neutrino background, in spite 
of the fact that we neglect contributions that should be included in a complete numerical study. 

While in the case that the neutrinos decay at rest, the photon flux vanishes when E > 6q, now it is possible to 
have a photon flux with energy higher than eo, albeit very small, since when a moving neutrino decays it can emit a 



Doppler blueshifted photon with energy higher than eg. From Eqs. (p4|), (35), and (ftOJ) we obtain 



cPFe I h Vi (t ) 1 rrii f E\ 2 \m; L E 



exp---l-, (47) 



dEdfl 4tt tH 3C(3) T V e o/ 2 T e ' 

valid for E ^> e , rrii 3> T , and tH not too small (so that the term exp— T/r does not dominate the fall of the 
spectrum) . 

To end this section we would like to point out that the photon spectrum obtained is not a thermal spectrum 
even if the neutrinos, that act as the photon source, are distributed thermally. The cosmic expansion and the decay 
kinematics render a photon distribution that is out of equilibrium. Later interaction of the decay products with 
matter could restore thermal equilibrium. However, we always assume that the decay happens after decoupling of 
matter and radiation, i. e., t 3> 300,000 yr, therefore the decay produced photons never attain thermal equilibrium. 



IV. CONCLUSIONS 



In this paper we have studied the contribution to the photon background of a cosmological density of decaying 
neutrinos. 

We first have adopted the assumption that neutrinos decay at rest and have investigated the fluxes in different 
cosmological scenarios: 1) a power-law cosmic expansion, 2) an open universe filled only with matter, and finally 3) 
a flat universe with matter and a nonzero cosmological constant A. We have calculated the exact formulae for these 
cases, and found the maximum of the spectrum as a function of the photon energy. 

The differences amongst the scenarios are not large. This result has the consequence that, when the observed 
photon spectrum is used to constrain the neutrino lifetime, the limit does not depend too much on the energy content 
of the universe. 

Second, we have studied the consequences for the photon background when the decaying neutrinos have a momentum 
distribution. Wc have worked out the relevant formulae for the scenarios 1), 2), and 3) mentioned earlier. We have 
also particularized our equations to the standard neutrino background. In all the cases, to compare with the photon 
spectrum calculated for neutrinos at rest, we have determined the first relativistic correction in our formulae. 

In the literature only the case of neutrinos decaying at rest in a critical matter dominated universe has been treated; 
this corresponds to Eq. ( |l4] ) . In this paper we have performed a general study of the photon spectrum that will allow 
to obtain reliable bounds on electromagnetic neutrino properties using photon background observations. Also, in 
the case of an eventual positive signal, the spectrum we have calculated would be helpful to obtain an observational 
determination of the lifetime r. 

A final comment is that we believe that our calculations may be useful when considering scenarios of radiatively 
decaying particles. For example, of the type considered by Sciama that might be the explanation for the observed 
ionization of intergalactic hydrogen. As another example, our results can be applied to the contribution to the 
photon background coming from a relic scalar or pseudoscalar particle density decaying into two photons of the type 
considered in jll],|l2|]. Of course one should introduce the appropriate changes in the equations of our paper. 
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APPENDIX A: DISTRIBUTION OF MOMENTA FOR A STABLE NEUTRINO COSMIC BACKGROUND 

In this appendix we shall study the momentum distribution of a relic cosmic background of stable neutrinos with 
mass m. We neglect the neutrino chemical potential setting it to zero although our results could easily be generalized 
to include a nonvanishing chemical potential. 

Neutrinos decouple from the primordial plasma of photons, electrons and baryons at td ~ 1 sec, when the temper- 
ature is Td ~ 1 MeV. Hereafter they freely stream without any further nongravitational interactions. Particle number 
conservation gives the following relation between the momentum distribution at time t > td and the momentum 
distribution at decoupling time t d 

f(k,t) = f(k d ,t d ), (Al) 

with the neutrino momenta related by k — kddd/ a, being a the cosmic expansion factor. Until time td the distribution 
function is given by the Fermi-Dirac thermal distribution, therefore 

f(k,t) = ) (A2) 

exp EdjTd + 1 



The energy at decoupling Ed — y m 2 + k 2 can be written as 




(A3) 

where the neutrino temperature T is related to the neutrino temperature at decoupling T — Td ad/ a. Writing 
everything together we obtain 




f(k,t)= l + ex P W- J + — . (A4) 



This distribution reduces to a thermal distribution only if the neutrinos are extremely relativistic, T > m, or 
nonrelativistic, T <C m, from decoupling time until present (for nonrelativistic neutrinos one has to redefine the 
temperature) The neutrino temperature T is related to the photon temperature by the relation T — (4/H) 1 / 3 ^. 
We are interested in a neutrino mass about m ~ 0.1 eV, hence the neutrinos are extremely relativistic at decoupling 
time but nonrelativistic at present because now T ~ 1.7 x 10~ 4 eV. Neglecting the term m/T d <C 1 in Eq. (|A3) we 
finally obtain 

f(k,t) = — !: . (A5) 

exp\k\/T(t) +1 

Note that, in general, this is not a thermal distribution since in the Fermi-Dirac distribution one has E = \/ m 2 + k 2 
which in our case has been replaced by \k\. Strictly speaking we cannot call the paramet er T neutrino temperature, 



but we use this term anyway since it does not cause any confusion. We make use of Eq. (A5) in Sec. Ill to calculate 
the spectrum of decay produced photons. 



APPENDIX B: CALCULATION OF THE FUNCTION T 



The integral T(y,z) defined in Eq. (37) can be solved in terms of elementary and special functions, 
consider separately the three cases introduced in Sec. ||. 



We shall 
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1. Power— law cosmic expansion 



When the Hubble function is given by the power-law (|Tl]) we have to calculate the following integral 

H T(y,z)= r ^ 1 (Bl) 
Ji+z v y/i + y v 

where n > 0. This integral can be related, by means of the change of integration variable u = y 2 v 2 /(l + y 2 v 2 ), to the 
incomplete beta function, which in turn is related to the hypergeometric function F(a,/3;7;it) = 2-Fi(a, P] 7; u) [ Of . 
For n not an even integer we obtain: 

H T{y, z) = - 

^V7T 



1 f Vl + y 2 (l + ^)M P / n l-n n y 2 (l + z) 2 



nl l + z / V 2 ' 2 ' 2'l + y 2 (l + z) 2 
For n an even integer, n > 0, the solution can always be given in terms of elementary functions 



F r(j/, z) = (-!)/ w , In -^— } + (i - z)n x 

\n n(n — 2) (n\\y ) 



(B3) 



Both Eqs. (|B2| ) and (B3) can be written in a more compact but less insightful way using the transformation formulae 
of the hypergeometric function [lg| ] 

H T(y,z) = — 1 =- — ?_ir(i,i ; !L±l; J -V (B4) 

Vy ' ^ n + 1 A /i + J/ 2( 1 + z )2 (l + z)» V 2 2 l + y 2 (l + 2) 2 7 ' k 7 

which holds for any n > 0. 

As particular examples of Eq. (|B2| ), we write down T(y,z) for a curvature and a flat matter dominated universe, 
n = 1, 3/2, respectively, 

H T(y, z )\ n=1 = -y H — , B5) 

l + z 



u T( y , z)| n=3/2 = -_l^r(i/4)V/ 2 

2(l + y 2 (l + z) 2 f 4 ( 3 1.1. y 2 (l + z) 



3 (l + z) 3 / 2 V 4 ' 4'4' l + j/ 2 (l + z) 2 J (B6) 

■J-r(i/4)V/ a + ^ 1 + 1/2(1 + z)a 

6Vf i / J ^ 3 (l + z)3/2 

-V 2 Ffarccosl^i±4, 1 V (B7) 
3^ I (l + y(l + z)) 2 'V2J 



The function F(cj>,k) is the incomplete elliptic integral of the first kind, see Sec. B2. In the limit y — > the above 
expressions become 

SoT(y, z)|„ =1 = flbt(«)|n=i - » + 0(y 2 ), (B8) 

H a T(y, z)\ n=m = H Q t(z)\ n=3/2 - ^=T(l/4) 2 y 3 / 2 + 0(y 2 ), (B9) 

where t(z)\ n is given by Eg. @ and r(l/4) 2 /60F = 1.2360 . . . 

As an example of Eq. (|B3|) one has a critical radiation dominated universe, n = 2, 
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^oT( y , z)| B=2 = —In tfi + z) + 2 (l + z) 2 • 



(BIO) 



In the limit y — > we find 



ffoT(y,z)|„ =2 = H t(z)\ n=2 + -y 2 \ny + 0(y 2 ), (Bll) 

with t(«)| n =2 given by Eq. @. 

As a final example of power law we consider a critical A dominated universe, n = 0. Directly integrating (Bl), or 
using (B4) (see representation of elementary functions in terms of a hypergeometric function in fl3|| ) we obtain 



H T(y,z)\ n=0 = In 



l + y/l + y 2 (l + z) 



y[l + z) 

When y —> this expression diverges because t(z) diverges for a critical A dominated universe. 



(B12) 



2. Open universe 



The integral to solve now is 



H T(y,z) 



dv 



l+z 



v y/fl k v 2 + n v 3 y/l + y 2 v 2 



This integral can be written as 



dv 



(B13) 



(B14) 



with P3 (v) a cubic polynomial, which is a particular case of elliptic integral JlJ] . By means of the variable change cn 



u = cn (u, k) = (v + a — A)/(v + a + A), whe re cn u is the Jacobian cosine amplitude u, a = Qk/^o, A 2 



1/y 2 , 



and k = (A + a)/2A, the elliptic integral (B13) can be written linearly in terms of elementary functions and the 
three fundamental incomplete elliptic integrals, of the first kind F, second kind E, and third kind II |l^ . [l5|l : 



y 1/2 ft . 



1 - 



V^o + W k y 2 



F(cf>,k) 



(n 2 + nly 2 ) 1/4 E(t,k) 



^n 2 + n 2 y 2 - n ky 



^lv 1/2 (ni + niy 2 ) 1 ' 4 
n 



-n(</),n 2 ,fc) 



An 



1 



3/2 



In 



n ky /i + y 2 (i + z) 2 + ^n k + n (i + z) 



/thy/i + y 2 (i + z) 2 - y/n k + n {i + z) 



(v^+W + n k y) y/n k + n (i + z) ^1 + 2,2(1 + ^2 



n k x /n 2 + n 2 y 2 + y(n k + n a (i + z)) 

We use the following definitions of the incomplete elliptic integrals |l4| ] 

d0 

F(cf>,k) = 
E{(f>,k) = 

n(0,n 2 ,/t) = 



l + z 



a/1 - k 2 



sin 



f d9 Vl - k 2 sin 2 6, 
Jo 



dO 



(1 - n 2 sin 2 6)\Jl - k 2 sin 2 9 



(B15) 

(B16) 
(B17) 
(B18) 
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The variables <j), k, and n of the three elliptic integrals are related to the physical variables of our problem: 

(B19) 



(n* + n (i + z)) y - ^n 2 + stiy 2 



n 2 



(r! fe + r» (i + z)) y + wTi 2 + n 2 y 2 

«*» Y (B20) 

1 /ng + n^fi+ " fct L = ) . (B2i) 



When y — > and Oq finite Eq. (B15) reduces to }L0 



iJ T(y, z) = H t(z) - J-r(l/4) 2 -Ly 3 / 2 + 0{y 2 ), (B22) 
6aA f^ /2 



where t(z) is given by Eq. (|l6|). 



For SIq — * and y finite Eq. (B15) gives 



H Q T(y, z) = H T(y, z)\ n=1 + ~ \H T(y, z)\ n=1 - In l+Vl+Z^+ fg ] fi + . . . (B23) 

2 V 2/(1 + *) / 

where the ellipsis means smaller contributions like fi 2 ,, f2 2 lnf2o, and so on. 

3. Flat universe with nonzero A 

In the last case that we study we have to cope with the integral 

H T(y,z)= / — - (B24) 

which is of the class 

dv 1 



1 + 2 



(B25) 



being P§(v) a five degree polynomial. Integrals of this sort can be solved as linear combinations of elementary functions 
and hyperelliptic integrals fll4f . Now there are five fundamental incomplete hyperelliptic integrals, two of the first 
kind, two of the second kind and one of the third kind. Since the solution in terms of the hyperelliptic integrals would 



not be particularly enlightening, we only study two limiting cases of (B24) 
For y — > and Qq finite we find 



H Q T(y, z) = H Q t(z) - -i-r(l/4) 2 -Ly 3/2 + 0(y 2 ), (B26) 
own oi /z 



ii 



where t(z) is given by Eq. |[|). In the regime Qq — > and y finite we obtain 



with 



H T(y, z) = H T(y, z)| n=0 - -Q l /o + . . . (B27) 

V 



F (- 1-v ^ 

6 = 1+ \ 3 ' 3 ' 3 ' 2 / = 1,199,., (B28) 

2 4 /3(l + V2) 2 /3 
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